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I.  INTRODUCTION 


In  recent  years,  the  Regenerative  Liquid  Propellant  Gun  (RLPG)  technology*  has  outpaced  the 
insight  into  the  complex  physics  of  its  underlying  injection/combustion  process.  A  major  conclusion 
from  earlier  research  using  relatively  simple  laboratory  test  fixtures^  was  that  a  more  advanced 
injection/combustion  facility^  had  to  be  constructed  if  further  insight  was  to  be  gained.  In  particular, 
an  emphasis  was  placed  on  constructing  an  injector  which  would  be  flexible  enough  to  allow  both 
fundamental  work  with  simple  jet  geometry,  e.g.  circular  (which  is  easier  to  model),  and  with  more 
practical  annular  jet  configurations  having  dimensions  approaching  those  of  common  30-mm  RLPG 
fixtures.  Thus,  requirements  for  the  injector  called  for  jets  with  thicknesses  of  1  to  3  mm  to  be 
injected  at  velocities  up  to  300  m/sec  into  ambient  pressures  exceeding  25  MPa.  A  desired  flexibility 
was  that  the  injection  velocity  could  be  varied  during  the  injection.  Clearly,  the  injector  to  be 
constructed  would  be  of  the  pressure  atomizing  type.  Common  laboratory  injectors  of  this  type  employ 
a  stepped  differential  area  piston.  The  larger  area  is  acted  upon  from  a  low  pressure,  large  external 
reservoir  while  the  smaller  area,  amplifying  the  reservoir  pressure,  acts  on  the  injected  liquid.  Such 
injectors  are  not  attractive  as  they  become  very  bulky  if  high  injection  velocities  are  to  be  sustained 
from  large  injection  ports,  such  as  annular  ones.  Also,  such  injectors  can  not  provide  varied  injection 
velocity  profiles.  From  a  safety  point  of  view,  when  dealing  with  liquid  propellant  (LP)  jet  combustion 
tests,  a  sudden  rise  in  ambient  pressure  would  result  in  a  drop  of  differential  injection  pressure,  hence 
risking  flash  back  into  the  injector. 

In  contrast,  an  RLPG  type  injector  is  a  compact  injector  with  high  mass  flow  rates  which 
inherently  responds  by  increasing  the  injection  pressure  when  the  ambient  pressure  suddenly  rises. 
Thus,  considering  the  alternatives  and  bearing  in  mind  that  the  type  of  injection  process  to  be  studied 
is  that  of  an  RLPG,  it  was  decided  to  construct  an  RLPG  type  injector  for  laboratory  use.  The  RLPG 
injector  utilizes  the  ambient  pressure  for  its  operation  but,  depending  on  its  design,  may  lack 
controllability  over  the  initiation  of  the  injection.  For  laboratory  use,  the  proposed  injector  had  to  be 
initiated  at  any  desired  initial  ambient  pressure.  This  was  solved  by  constructing  an  injector  that 
employs  a  pressurized  gas  chamber  as  a  counter  balance  against  the  ambient  pressure.  Control  over 
injection  initiation  and  injection  velocity  profile  is  achieved  by  controlling  the  pressure  discharge  from 
the  gas  chamber.  This  report  details  the  structure  and  performance  of  the  injector  which  was 
constructed  for  use  in  the  BRL  Spray  Research  Facility^  (SRF).  The  injector  features  compactness, 
rapid  response,  circular  as  well  as  annular  jet  configuration,  good  controllability,  and  safety  of 
operation.  To  date,  the  injector  has  been  tested  successfully  in  all  configurations  at  25  MPa  ambient 
pressure.  LP  combustion  tests-*  have  been  conducted  featuring  volumetric  flow  rates  of  up  to  1000 

cc/sec  for  total  injected  volumes  up  to  92  cc. 

The  purpose  of  this  report  is  to  serve  as  a  manual  for  the  operation  of  the  injector  and  as  a 
reference  in  future  reported  work  from  the  SRF.  The  reader  is  also  given  the  computer  code  which  is 
used  for  the  prediction/validation  of  the  injector  performance.  The  information  contained  herein  may 
enable  researchers  to  construct  similar  injectors  and  adapt  them  to  their  own  uses.  Presently,  a  patent 
is  pending  on  the  injector. 
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n.  INJECTOR  STRUCTURE  AND  OPERATION 


A  simplified  assembly  drawing  (drawn  to  scale)  of  the  injector  is  shown  in  Figures  la,  lb,  and 
lc.  An  exploded  view  of  the  injector  is  depicted  in  Figure  2.  The  major  parts  and  features  of  the 
injector  are  identified  on  both  Figures  1  and  2.  The  injector  has  three  modes  of  operation  K,  J,  and  L. 

The  basic  mode,  K,  generates  a  circular  jet  while  modes  J  and  L  generate  annular  jets.  All  modes 
consist  of  four  major  part  assemblies  l,  2,  3,  4  and  chambers  ’A,  B,  C’  (and  ’D’  in  case  of  modes  of  J 
and  L)  which  contain  the  fluids  involved  with  the  operation  of  the  injector.  Chambers  ’A’  and  ’D’ 
contain  liquid  while  chambers  ’B’  and  ’C’  contain  gas  (typically  nitrogen).  Chamber  ’B’  is  shown  in  the 
closed  position  in  Fig.  la.  All  other  parts  are  assembled  into  parts  1,  2,  3,  and  4.  Part  1  is  the  outer 
housing  of  the  injector  and  is  sealed  against  the  bottom  of  the  test  chamber  walls  30,  thus  forming  the 

base  of  test  chamber  *T’  into  which  the  liquid  is  injected  through  port  ’j’.  Part  2  forms  the  base  of  the 

,  .  * 

injector  into  which  all  gas  supplies  are  connected.  Part  3  (or  3  in  modes  J,  L)  is  the  injector  outer 
* 

piston.  Part  4  (or  4  in  mode  L)  is  the  injector  inner  piston.  Parts  1,  2,  and  4  are  stationary  while 
part  3  is  mobile.  The  liquid  to  be  injected  is  contained  in  chamber  ’A’  which  forms  between  3  and  4. 
The  maximum  injection  volume  is  92  cc  for  mode  K  and  45  cc  for  modes  J  and  L.  Piston  4  is  held  on 
top  of  2  by  the  pressure  P^  in  ’A’.  The  injector  parts  which  are  exposed  to  the  liquid  are  made  of 
stainless  steel  316.  Other  parts  are  made  of  stainless  steel  17-4PH.  The  injector  is  designed  for  70 
MPa  differential  injection  pressure.  The  external  gas  controls  for  the  injector  will  be  detailed  in 
reference  5. 

The  principle  of  operation  of  the  injector  is  similar  to  that  of  of  the  RLPG.  For  its  operation, 
it  is  required  that  a  pressure  P^-.  be  established  in  ’T\  Due  to  the  differential  areas  of  part  3,  a 
pressure  P^  is  formed  in  ’A’  which  causes  the  liquid  in  ’A’  to  be  injected  through  ’j’.  With  the  absence 
of  pressures  in  ’B’  and  ’C’,  the  differential  area  of  3  results  in  P^  =  2.29  Pc  in  Modes  K  and  J,  and  in 
PL  =  2.00  P^  in  Mode  L.  Unlike  the  RLPG,  the  subject  injector  has  two  gas  chambers  ’B’  and  ’C’ 
which  serve  respectively  to  augment  or  arrest  the  motion  of  3  by  changing  the  force  balance  on  3. 
During  injection,  the  piston  3  moves  downward  which  decreases  the  volumes  in  ’A’  and  ’C’  (and  ’D’  in 
mode  L)  and  increases  the  volumes  in  ’B’  and  ’T\  Adding  pressure  in  ’B’  and  T’  while  decreasing  the 
pressure  in  ’C’  will  result  in  an  increase  in  P^  and  higher  injection  velocity.  The  detailed  physical 
formulation  of  the  injector  operation  is  given  in  Section  HI.  Following  are  the  details  concerning  the 
injector  subassemblies  and  operation.  Mode  K  is  discussed  first. 

Housing  1  contains  the  exhaust  port  for  ’T’  -  ’exh\  The  pressure  in  ’C’  is  monitored  by  a 
601B1  Kistler  pressure  transducer  -  ’p’.  Three  pins  5  provide  a  stop  for  3  when  30  is  lifted  away  from 
1.  When  assembled  into  ’T’,  3  is  stopped  against  30  by  6.  Pins  5  also  prevent  3  from  rotating  while  in 
longitudinal  motion  as  they  protrude  into  vertical  tracks  machined  into  3.  Part  2  contains  all  the  gas 
ports  leading  into  chambers  ’B’  and  ’C\  ’C’  is  filled  with  gas  via  the  port  ’chg’  and  discharged  via  port 
’deg’.  A  port  ’aug’  is  used  to  fill  chamber  ’B’  (and  also  to  discharge  it  after  injection).  Another  port, 
not  shown,  is  used  to  relieve  gas  or  liquid  leaks  which  may  occur  across  the  seals  between  4  and  10  and 
between  2  and  4. 

Three  subassemblies  are  mounted  into  2.  The  first  of  which,  15,  is  used  to  record  the  piston 
motion.  A  tapered  shaft  14  screws  into  3  and  passes  through  2  and  15.  Camshaft  mechanism  32 
changes  proximity  to  a  Hall  effect  magnetic  sensor  16  (Wolff  Microjector)  which  results  in  a  voltage 
output  polynomially  fitted  to  the  vertical  position  of  piston  3.  Thus,  piston  velocity,  chamber  volumes, 
and  therefore  injection  velocity  can  be  easily  calculated  as  a  function  of  time. 


2 


I 

Fig.  2.  Exploded  view  of  injector  parts. 
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The  second  subassembly  is  an  ultra-fast  exhaust  valve  connected  to  port  ’deg*.  Part  19  is 
connected  to  port  ’deg’  by  the  stainless  steel  tube  shown  in  Fig.  2.  The  valve  has  a  differential  area 
stepped  piston  20,  an  upper  cup  19,  housing  18  and  a  bottom  plug  21.  In  its  closed  position,  piston  20 
is  pushed  into  19  such  that  21  seals  against  18.  A  low  pressure  (typically  0.8  MPa)  is  introduced  from 
the  laboratory  air  supply  through  ’g’  into  G.  The  piston  20  then  presses  21  onto  18  thus  strengthening 
the  seal  against  18.  The  pressure  in  G  is  sufficient  to  overcome  the  pressure  Pjj  (the  pressure  in  *C’) 
which  exists  in  19  and  exerts  on  20.  Discharge  of  the  pressure  in  G  through  ’g’  results  in  20  moving 
under  PH  and  instantly  unseating  21  from  18  which  results  in  a  prompt  discharge  of  the  gas  in  G 
through  the  gap  forming  between  21  and  18.  As  20  accelerates  downward,  its  upper  step  departs  19 
and  the  gas  from  *C*  is  suddenly  discharged  through  port  ’e’.  Various  orifices  can  be  mounted  in  19  in 
order  to  control  the  discharge  mass  flow  rate  from  ’C’  through  19.  This  controls  indirectly  the 
injection  velocity  by  virtue  of  regulating  PH-  In  its  fastest  release  mode  (i.e.  no  restricting  orifices), 
the  valve  provides  an  instant  (sub-millisecond)  opening  of  .200  inch  effective  diameter.  Such  an 
ultra-fast  valve  is  necessary  in  order  to  provide  the  required  wide  range  of  injection  flow  rates  and 
velocity  profiles. 

The  third  subassembly  mounted  into  2  serves  a  dual  role:  sealing  the  injection  port  ’j*  prior  to 
injection  and  providing  a  means  for  measuring  both  port  ’j’  opening  and  pressure  P, .  (A  Kistler 
607C3  pressure  transducer  17  is  mounted  into  plug  13  and  measures  PL.)  A  grease  column  ’GR’ 
between  17  and  piston  12  transmits  to  17  the  pressure  which  is  exerted  on  12.  Prior  to  injection,  the 
piston  3  is  in  its  uppermost  position  as  defined  by  6  exerting  on  30  (a  spacer,  not  shown,  may  be 
placed  between  6  and  30  which  results  in  a  smaller  initial  ’A’  and  ’C’  and  larger  ’B’).  The  gas  in  ’C’ 
enters  the  void  between  12  and  11  via  passage  ’b’.  Pilot  shaft  10  lifts  under  PH  until  it  seals  against  8 
thus  preventing  injection.  The  area  ratio  across  10  is  such  that  a  P^  >  2.25  Pjj  is  required  in  order  to 
unseat  10  from  8.  When  the  pressure  Pjj  is  discharged  through  ’deg’,  10  retracts  from  the  injection 
port  and  exerts  force  on  12  via  11  which  translates  into  pressure  measured  by  17.  Therefore  prior  to 
injector  port  opening,  17  measures  PH'  As  11  contacts  12,  17  measures  P^,  though  indirectly.  When 
11  contacts  12,  Pjj  is  exerted  equally  on  the  bottom  of  10  and  on  the  top  of  12,  thus  the  net  force  felt 
by  17  is  solely  due  to  P^  exerting  on  the  portion  of  10  above  its  seal  in  ’A’.  Prior  to  injection,  10 
serves  to  insulate  the  combustible  liquid  in  ’A’  from  the  high  temperature  gas  in  *T\  It  also  delays 
piston  3  motion  until  the  pressure  PH  drops  below  P^  /  2.25,  thus  resulting  in  higher  initial  injection 
velocity.  Part  10  exposes  ’j’  in  a  controlled  manner,  thus  resulting  in  a  smooth  injection  initiation. 
The  injector  can  also  be  operated  with  shaft  10  (Fig.  lc)  instead  of  10.  This  will  be  referred  to  as 
"no  shaft"  opening  mode  (as  opposed  to  "pilot  shaft”).  Piston  3  will,  in  this  case,  start  its  motion  when 
PH  <  PC  /  1.71. 

In  order  to  refill  chamber  ’A’  with  liquid,  30  is  lifted  from  1,  and  3  is  pushed  to  its  desired 
position  by  applying  some  pressure  in  *C\  Cap  6  and  upper  plug  7  are  removed  as  well  as  teflon 
barrier  9.  An  injector  plug  8  (with  a  selected  diameter  for  injector  prat  ’j’)  is  screwed  onto  4.  The 
liquid  is  poured  into  volume  ’A’  which  includes  passage  ’a’  and  the  space  around  10.  A  bleed  hole  (not 
shown)  exists  in  7.  Parts  6,  7  and  9  are  reassembled  and  the  bleed  hole  plugged.  During  injection,  3 
moves  downward,  thus  displacing  the  liquid  from  ’A’  via  ’a’  around  10  and  injecting  it  into  T*  via  ’j’. 
Barriers  9  serve  as  piston  guides  and  also  as  insulation  between  the  liquid  and  moving  seals.  This  is  a 
safety  precaution  for  the  cases  when  the  liquid  is  LP.  For  spray  combustion  tests  a  water  layer  may  be 
placed  on  top  of  the  LP  in  ’A’  below  the  upper  barrier  9.  A  thin  mylar  barrier  may  be  used  to 
separate  the  LP  and  water.  As  the  LP  is  displaced  from  ’A’,  the  water  flushes  ’A’  and  prevents  ignition 
in  ’A’  when  the  piston  3  is  brought  to  a  halt  by  rubber/teflon  shock  absorber  31. 
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Mode  J  (Figure  lb)  operates  virtually  the  same  as  mode  K.  The  piston  3*  is  similar  to  3  in 
mode  K,  it  is  just  longer.  The  injection  port  ’j’,  formed  between  24  and  7  ,  is  annular  and  unlike 
mode  K,  it  is  not  sealed  from  *T’  prior  to  injection.  Piston  3  is  held  against  30  by  split  ring  22.  A 

floating  piston  23  separates  water  in  ’D’  horn  the  injected  liquid  in  ’A’.  It  also  supports  the  motion  of 

3  along  4.  The  annulus  thickness  is  varied  by  replacing  7  with  a  part  of  a  different  inner  diameter. 
The  operation  of  mode  J  results  in  an  injection  pattern  much  like  RLPG  concept  VI  where  a  center 
body  24  is  exposed  as  the  piston  3  retreats.  Piston  3  will  start  its  motion  when  PH  <  pC  /  l  ^- 

Mode  L  (Figure  lc)  was  constructed  in  order  to  simulate  RLPG  concept  VI-C.  The  injection 
port  is  also  annular.  Here  no  center  body  is  exposed  as  the  piston  retreats.  Cap  26  has  a  geometry 
similar  to  the  one  found  in  a  30  mm  RLPG  concept  VI-C.  Cap  26  is  mounted  on  pedestal  25  which  is 
held  into  3  by  means  of  snap  rings  27  and  spacer  29.  Water  fills  volume  ’D’  while  the  injected  liquid 
fills  volume  ’A’.  The  motion  of  3  displaces  the  water  from  ’D’  which  pushes  23  upward,  thus 
injecting  the  liquid  from  ’A’.  Plugs  28  delay  injection  as  they  provide  an  impediment  to  the  flow  of 
the  water  from  ’D’  through  the  holes  in  the  bottom  of  25  onto  the  bottom  of  23.  Thus,  most  of  Pjj  can 
be  removed  before  a  significant  amount  of  water  is  displaced  from  ’D’  through  25.  The  geometry  of  ’j’ 
can  be  varied  by  using  different  thickness  spacer  29.  Also  plugs  24  and  1*  from  mode  J  can  be  placed 
in  lieu  of  26  and  7  thus  resulting  in  a  simpler  geometry  as  in  mode  J.  Mode  L  can  be  operated 

without  the  plugs  28.  With  the  plugs,  piston  3  will  start  its  motion  when  Pjj  drops  to  some  value 

between  Pc  /  1.71  and  Pc  /  3.77. 

There  are  various  ways  for  controlling  the  injection  process.  They  are  discussed  in  Section  IV 
together  with  examples  of  the  predicted  performance  of  the  various  modes  of  operation. 
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HI.  FORMULATION  OF  INJECTOR  OPERATION 


In  this  section,  the  operation  of  the  injector  in  its  various  modes  of  operation  is 
formulated.  The  following  assumptions  are  made: 

a)  the  liquid  is  considered  incompressible. 

b)  the  gas  flow  is  isentropic. 

c)  the  friction  forces  are  negligible. 

d)  the  liquid  flow  is  considered  steady,  non-viscous  and  one  dimensional. 

Assumption  ’a’  simplifies  the  formulation  a  great  deal  and  is  justified  considering  the  type  of 
liquids  and  pressure  range  under  consideration.  Assumption  ’b*  introduces  small  errors  in  long 
injection  processes.  These  errors  can  only  be  determined  experimentally.  Assumption  *c*  is 
better  for  the  higher  pressures;  its  validity  has  to  be  determined  experimentally.  Assumption ’d’ 
simplifies  the  solution  and  is  justified  based  on  the  operating  conditions  considered.  An 
estimated  discharge  coefficient  is  introduced  in  order  to  account  for  assumption  ’d\  The 
uncertainties  in  Cjj  mask  those  of  assumptions  ’a’,  *b\  and  ’c*.  Reasonable  deviations  from  the 
above  assumptions  were  simulated  in  the  equations,  and  the  solutions  indicated  less  than  5% 
deviation  in  injector  performance.  Therefore,  for  the  purpose  of  predicting  the  injector 
performance,  the  above  assumptions  are  justifiable. 

* 

A.  Equations: 

Balance  of  forces  on  the  piston: 

PAUAWJ+  PCA0  =  PlAI  +  PHAH+  MP'd't  ^ 


Mass  continuity  of  the  liquid  during  injection: 


UpA.=  UjAj 
Up(A,  +  At)  =  UjAt 


for  modes  K,  } 
for  mode  L 


Steady  Bernoulli  Equations  for  the  liquid  flow: 

J*fl2 


P  -  P  = 

rL  rC 


C  2 


(2) 


(3) 


Isentropic  gas  equation  of  state  in  gas  chambers: 

p.  =  c.pj 


C. 


y-i 

_  MP..  ’ 

~*rxr 


(4a) 

(4b) 


irtiere  .  stands  for  H  ^  c  (corresponding  to  gas  chambers  C,  B.  and  T  In  Fig.  I) 
and  stands  for  jjq  5  respectively. 


*Symbols  defined  in  Nomenclature,  page  23. 
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(5) 


For  conciseness,  the  exact  forms  of  Eqs.(8)  are  not  given  here.  Their  forms  can  be  determined 
from  the  Program  INJFLOW  given  in  Appendix  A. 
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The  initial  conditions  for  Eqs.(8)  are: 


PH(t=0)  =  pm  (t=0)  =  Pffl0  pc  (t=0)  =  ps  (9) 

VH  (t=0)  =  Vm  ^  (t=0)  =  0 


Notes: 

1)  Equation  (8b)  holcj  for  the  case  when  the  test  chamber  is  interconnected  to  a  large  capacity 
pressure  reservoir  Pg. 

2)  A  different  mode  of  operation  results  when  is  prescribed  in  a  given  manner.  An 
exponential  pressure  rise  is  considered  (simulating  a  gun  cycle  after  ignition,  or  injection  during 

a  rapid  pressurization  of  the  test  chamber  to  a  pressure  P§)  in  the  form: 

Pc  =  Ps(l  -  eaexP‘)  (10) 

In  this  case  Eq.  (9)  replaces  Eq.  (8b)  and  Eqs.  (8)  are  reduced  to  4  ODE’s  . 


Condition  fa*  injector  initiation  (injector  opening)  will  now  be  described: 


Note  that  Eq.  (1)  is  valid  only  when  the  hold  pressure  has  vented  to  a  point  when  the 
test  chamber  no  longer  exerts  a  restraining  face  on  the  face  of  the  piston.  It  is  only  then  that 
the  piston  may  start  its  motion,  thus  initiating  the  injection.  The  piston  however  may  be 
blocked  from  movement  due  to  excess  liquid  pressure  build  up  behind  the  plugs  (28  in  Figure 
lc)  in  Mode  L  or  behind  the  pilot  shaft  (10  in  Figure  la)  in  Mode  K.  As  the  hold  pressure  is 
vented  further,  the  plugs  and/or  the  pilot  shaft  recede  thus  triggering  the  injection.  In  order  to 
account  for  the  various  modes  of  operation,  Eqs.  (8)  have  to  be  solved  with  respect  to  the 
various  modes  of  injector  opening  (ie.,injector  initiation  at  t=tQP).  Therefore  Eqs.  (8)  are 
subjected  to  the  following  constraints: 


-ajffcctop)  » o 


(lla) 


PL  (t<top)  -  Pc 
p  ,  ^aiAuj  +  pc  ”  AP 

HOP  s  JTh 

p  „  pcAo  +  pauaau 
hop  -  z.z5a{  +  ah 


PH  "  PH0P  ,(hen  1  '  *0P 

except  vhen  opening  mode  of  pilot  shaft  or  plugs 


for  Node  K  vlth  pilot  shaft 


(lib) 

(11c) 

(lid) 


HOP 


2  P r 


Agfro  +  AT>  ~  (A,  +  AT  - 


ApS)  (Aj  +  at  ~ahs) 


amsah 


for  Mode  L  vlth  plugs  (lie) 


B.  Solution: 


A  fourth  order  Runge-Kutta  method  is  used  to  solve  Eqs.  (8).  The  equations  are  stiff 
ODE’s  and  extensive  logic  is  used  fa  controlling  the  time  steps  in  order  to  minimize  their 
number  and  yet  assuring  solution  convergence  and  accuracy.  The  ODE’s  are  made  less  stiff  if 
the  piston  inertia  term  in  Eq.  (1)  is  neglected  which  renders  Eq.  (8e)  algebraic.  Therefore,  the 
value  of  the  inertia  term  is  constantly  checked  during  the  solution  and  it  is  neglected  if  its 
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magnitude  is  much  smaller  with  respect  to  the  other  terms  in  Eq.  (1)  (or  their  equivalent  in  Eq. 
(8d)).  If  the  inertia  term  is  neglected  then  only  4  ODE’s  are  solved  (or  3  ODE’s  if  Eq.  (10) 
applies).  As  Eq.  (11)  indicates,  only  1  ODE  is  solved  for  t  <  tQp.  For  the  simulation  of  the 
injector  operation,  a  Microway  FORTRAN  program  INJFLOW  was  written  (APPENDIX  A)  and 
run  on  a  386  based  PC-AT. 

Input: 

The  various  modes  of  operation  of  the  injector  and  its  operating  parameters  are 
interactively  input  at  the  terminal  according  to  the  flow  chart  given  in  APPENDIX  B. 

Output: 

The  various  injector  performance  variables  are  displayed  on  the  monitor  in  a  tabular 
form  during  die  solution  execution.  They  are  also  simultaneously  stored  in  ASCII  files  on  the 
system  disk.  The  files  automatically  assume  coded  names  in  accordance  with  the  input  for  the 
runs.  All  file  names  are  eight  characters  long,  plus  a  three  digit  suffix.  The  first  character 
declares  the  type  of  injection  being  modeled.  It  is  either  ’K’  for  a  circular  injection  (mode  K) 
or  ’L’  or  T  for  an  annular  injection  (mode  L  or  J  respectively).  The  second  character  is  either 
’P’  for  the  preset  test  chamber  pressure  mode  (Eq.  9),  or  ’E’  for  the  exponential  test  chamber 
pressure  mode  (Eq.  10).  The  third  through  fifth  characters  are  either  the  hold  discharge  orifice 
diameter  in  thousandths  of  an  inch,  or  are  ’HHH’  if  constant  hold  pressure  is  kept.  The  sixth 
character  has  different  meanings  depending  on  whether  the  injection  was  circular  or  annular.  If 
circular,  it  is  either  ’S’  for  pilot  shaft  (10)  present,  *N*  for  no  pilot  shaft  present,  or  ’A’  for  an 
arbitrary  opening  hold  pressure  (a  shaft  is  present).  If  annular,  it  is  ’P’  for  using  plugs  (28),  ’N’ 
for  not  using  plugs,  or  ’A’  for  an  arbitrary  opening  hold  pressure  (plugs  are  present).  The 
seventh  character  gives  the  value  for  CD  in  tenths,  and  reads  ’O’  if  CD  »  1.  The  eighth 
character  gives  the  liquid  injected.  It  is  ’P’  for  propellant  (LP),  ’W’  for  water,  or  ’M’  few  a 
miscellaneous  liquid.  The  three  digits  in  the  suffix  give  the  injection  orifice  diameter  in 
thousandths  of  an  inch  for  circular  injections,  or  the  injection  orifice  diameter  minus  an  inch 
for  annular  injections  (outer  diameter  for  mode  J,  inner  diameter  for  mode  L). 
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IV.  INJECTOR  PERFORMANCE 


In  this  section,  the  predicted  performance  of  the  injector  is  given  for  selected  operating 
conditions.  It  is  not  intended  here  to  cover  the  entire  performance  envelope  of  the  injector,  but 
rather  to  demonstrate  the  trends  of  the  performance  as  calculated  based  on  the  formulation  of 
Section  m.  This  will  cue  the  user  on  the  feasibility  of  desired  injection  velocity  profiles.  The 
theoretical  predictions  are  given  in  Figures  3  to  8.  Actual  test  results  and  comparison  to  theory 
are  given  in  Figure  9.  In  all  Figures,  the  file  name  corresponding  to  the  operational  parameters 
used  are  given  at  the  top  (e.g.,  KP200N9P.020  in  Fig.3);  thus,  the  parameters  can  be  deciphered 
as  explained  in  Section  HI. 

A.  Predicted  Performance 

In  Figs.  3-8,  the  hold  volume  is  set  such  that  the  volume  injected  is  approximately  43 
cc.  Also  the  initial  hold  pressures  are  set  such  that  they  are  only  slightly  above  the  minimum 
required  to  prevent  injection.  The  augment  pressure  is  always  applied  only  after  injection 
initiation.  In  all  cases  LGP1846  is  the  liquid  injected.  (Lighter  liquids  will  result  in  higher 

injection  velocities.) 

The  injector  performances  for  mode  K  (circular  orifice)  are  given  in  Figs.  3-5  for  the 
range  of  ambient  pressures  which  will  be  tested  in  the  SRF.  In  each  Figure,  injection  velocities 
were  calculated  for  cases  with  and  without  augment  pressure.  (Note  that  the  augment  pressure 
should  always  be  lower  than  the  ambient  pressure;  otherwise  Part  4  will  unseat  from  Part  2. 
For  identification  of  Parts  and  Volumes,  see  Fig.  1.)  Figures  3  and  4  cover  extremes  in  the 
injector  orifice  size  (small  and  large,  respectively)  when  the  hold  discharge  orifice,  A^,  is  the 
largest  possible  size.  A  large  AyH  results  in  faster  rates  of  PH  decrease  and  P^  increase.  For 
the  small  injection  orifice  (Fig.  3),  the  injector  piston  (Part  3)  moves  slowly  enough  such  that 
PH  drops  to  atmospheric  value  very  early  in  the  injection  period  which  results  in  steady,  high 
injection  velocities.  Also,  the  augment  volume  (’B’)  expands  at  a  slow  rate  such  that  PAy  is 
sustained  in  its  initial  set  value  despite  the  small  inlet  area  AyAjj.  In  contrast,  for  the  larger 
injection  orifice  (Fig.  4),  the  injection  piston  moves  faster  which  recompresses  the  gas  in  ’C* 
such  that  Pj|  never  drops  to  atmospheric  value  during  the  injection  period.  This  results  in  an 
accelerating  injection  velocity  during  injection  and  the  top  injection  velocity  is  smaller  than  the 
one  achieved  with  the  smaller  injector  orifice  for  the  same  operating  pressures.  Notably,  PAU 
drops  as  the  gas  supply  through  Aj^j  cannot  keep  up  with  the  rapidly  expanding  ’B\  As  can 
be  seen  from  Fig.  4,  with  the  proper  level  of  PAy,  the  injection  velocity  profile  can  be  made 
flatter.  In  Fig.  5,  a  medium  sized  injection  orifice  is  used  with  a  small  Ayjj.  This  results  in  a 
long  duration,  steady  albeit  low  injection  velocity.  The  velocity  stabilizes  early  as  the  volume 
’C’  decreases  at  a  rate  which  compensates  exactly  for  the  discharge  of  gas  from  it  through  the 
small  Ayjj.  When  augment  pressure  is  used,  the  initial  piston  velocity  is  high  but  settles  down 
as  both  PAU  and  PH  stabilize. 

Note  that  Figs.  3-5  deal  with  injector  opening  mode  of  no  pilot  shaft  (Part  10).  As  seen 
in  the  detailed  inset  of  Fig.  S,  when  PAy  is  applied,  P^  is  compressed  to  values  above  its 
initial  value  which  would  cause  Part  10  in  a  pilot  shaft  mode  to  seal  the  injection  orifice  and 
thus  terminate  the  injection  prematurely.  In  Figs.  4  and  5,  the  inset  details  are  for  cases  when 
PH0  ‘s  561  weU  above  the  value  required  to  prevent  injection.  As  P^q  is  discharged,  injection 
starts  and  is  then  boosted  further  by  the  application  of  PAU'  In  Fig.  4,  the  injection  velocity 
profiles  could  have  been  made  flatter  had  a  pilot  shaft  opening  mode  been  used.  This  is 
because  in  this  mode  the  injection  velocity  starts  at  a  value  above  100  m/s. 
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The  cases  covered  in  cases  Figs.  3-5  are  examples  which  show  the  trends  in  expected 
performance  with  respect  to  variables  such  as  the  injection  orifice  size,  pressures  on  the  injector 
piston,  and  the  hold  discharge  orifice.  As  shown,  for  small  injection  orifices,  the  injection 
duration  is  long  and  the  hold  volume  ’C’  changes  slowly.  For  such  a  case  it  is  possible  to  hold 
Ppj  constant  at  a  preset  value  by  means  of  a  regulator.  Injection  may  be  initiated  by  the 
application  of  PAU  which  is  the  case  shown  in  Fig.  6.  Injection  can  also  be  initiated  by 
increasing  the  ambient  pressure  as  shown  in  Fig.  7  for  an  exponential  ambient  pressure  rise. 

In  Fig.  7,  P^y  is  atmospheric  and  only  a  small  PH0  is  required.  Obviously,  variations  on  these 
themes  may  be  used  in  order  to  achieve  any  desired  injection  profile. 

As  Figs.  3-8  indicate,  it  is  reasonable  to  expect  injection  velocities  from  200  m/s  to  over 
300  m/s  for  circular  jets  of  LGP1846  which  range  in  diameter  from  3.2  mm  (0.125  inch)  to 
below  0.5  mm  (0.020  inch)  respectively.  For  the  thickest  jets,  the  top  injection  velocity  is 
achieved  only  toward  the  end  of  the  injection  period.  The  injection  velocities  attainable  with 
the  annular  jets  are  lower  than  with  the  circular  jets.  This  is  because  the  annular  jets  have 
larger  cross  sectional  areas  which  require  higher  injector  piston  velocities  in  order  to  achieve  the 
same  high  injection  velocities.  This  in  turn  recompresses  the  gas  in  the  hold  volume,  thus 
impeding  injection. 

In  mode  L,  shown  in  Fig.  8,  the  opening  injector  pressure  PjjQp  can  be  controlled  by 
the  use  of  plugs  (Part  28).  For  example,  for  Py  =  27.7  MPa,  without  the  u.  •  of  plugs,  a  top 
velocity  of  25  m/s  is  obtained  with  an  annular  jet  1  mm  (0.040  inch)  thick  and  25.4  mm  (1.0 
inch)  ID.  For  the  same  jet.  with  the  use  of  the  plugs,  injection  velocity  starts  at  150  m/s  and 
drops  to  90  m/s  during  the  5  ms  injection  period.  A  set  of  newly  designed  plugs  are  expected 
to  open  at  Pjjqp  =  l-4^  MPa  which  will  result  in  190  to  175  m/s,  over  a  3  ms  injection.  Note, 
that  Py  drops  because  the  ambient  gas  is  not  provided  quickly  enough  to  allow  for  the 
expansion  of  ’A’  as  the  injector  piston  withdraws  from  the  bottom  of  ’A*.  The  performance  of 
Mode  J  injection  which  lacks  the  plugs  is  rather  similar  to  the  performance  of  Mode  L  without 
the  plugs.  The  maximum  injection  velocity  is  limited  by  the  minimum  Pjjq  which  is  required 
to  hold  the  piston  against  Py,  and  by  the  maximum  available  discharge  area  Ajjj.  If  the 
injector  is  operated  as  a  true  RLPG  injector,  i.e.,  with  minimal  P^q  under  exponentially 
increasing  Py,  then  injection  velocities  in  the  range  of  300  m/s  are  achievable  even  with  3  mm 
(0.120  inch)  thick  -  25.4  mm  (1.0  inch)  ID  annular  jets.  Obviously,  the  injection  durations  will 
be  very  short  which  means  that  Py  has  to  be  applied  very  rapidly.  This  can  be  achieved  only 
by  employing  a  high  output  gas  generator  in  lieu  of  the  present  nitrogen  gas  supply.  The 
injector  performance  for  thick  jets  can  be  improved  significantly  if  it  is  operated  with  helium 
in  lieu  of  nitrogen  but  the  operation  will  be  less  practical. 

B,  Experimental  Performance 

The  predicted  performance  of  the  injector  was  validated  in  a  series  of  tests.  The  test 
shown  in  Fig.  9  is  typical  for  mode  K  (circular)  injection  with  the  use  of  the  pilot  shaft  (Part 
10).  P£  was  preset  at  28  MPa  and  P^q  was  set  at  14.5  MPa  (time  A).  Experimentally,  the 
pilot  shaft  retracted  from  the  injection  port  at  about  9.5  MPa  (time  B)  thus  opening  the  injector 
and  initiating  the  injection.  As  explained  in  Section  2,  from  time  A  to  B,  transducer  17 
measures  Pjj  directly  and  thus  its  output  is  the  same  as  transducer  ’p’  (Fig.  1).  During  the 
injection  (time  B  to  C),  ’p’  continues  the  measurement  of  PH  while  transducer  17  measures  P^ 
via  shaft  10.  The  output  of  17  during  injection  has  to  be  multiplied  by  a  2.25  area  factor,  hence 
the  new  scale  for  P^  at  time  B,  The  experimental  curve  was  calculated  from  the  magnetic 
sensor  (Part  16)  data.  (About  78  cc  of  water  was  injected.)  P^  dropped  significantly  during 
the  injection  due  to  the  expansion  of  ’A’.  (The  gas  supply  level  was  too  low  in  the  experiment.) 
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KPHHHN9P.020 


Fig.  6  Initiation  of  injection  by  application  ol  augment  pressure  while 
maintaining  the  hold  pressure  at  a  preset  value. 


Fig.  7  Initiation  of  injection  by  pressurization  of  the  test  chamber  while 
maintaining  the  hold  pressure. 
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75A0W.080 


Comparison  between  theoretical  analysis  and  experimental  data  on  the 
operation  of  the  injector. 


Experimentally,  the  most  precise  data  are  for  PH.  PC.  and  the  times  A,  B,  and  C.  The  least 
precise  data  are  for  P^.  For  validation  of  the  theoretical  model,  the  code  INJFLOW  was  run 
with  the  experimental  opening  pressure  of  the  injector  (Prop  01  t"ne  as  an  input  As 
evidenced  from  Fig.  9,  the  agreement  between  the  experiment  and  theory  is  quite  good.  As 
expected,  the  model  predicted  better  performance  when  a  discharge  coefficient  of  one  was  used 
and  no  losses  were  taken  into  account 
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V.  CONCLUDING  REMARKS 


1.  Based  on  model  predictions  and  supported  by  experiments,  the  injector  performance  is 
adequate  for  the  studies  of  RLPG  type  injection  processes.  With  circular  geometries,  the  full 
range  of  RLPG  jet  velocities  can  be  researched  with  mass  flow  rates  up  to  3  kg/s.  With  annular 
geometries  approximating  a  30- mm  RLPG,  mass  flow  rates  up  to  10  kg/s  may  be  obtained  but 
with  attendant  injection  velocities  below  200  m/s.  With  circular  jets  up  to  92  cc  of  liquid  can 
be  injected,  while  with  annular  jets  injection  volume  is  limited  to  45  cc. 

2.  Operation  with  the  thinner  (2  mm  (0.080  inch)  in  diameter)  circular  jets  offers  long  injection 
durations  with  steady  injection  velocities  which  are  needed  for  statistical  analysis  of  particle 
sizing  data,  generated  by  Phase  Doppler  Particle  Analyzer  measurements  of  the  spray. 

3.  The  use  of  the  pilot  shaft  opening  mode  is  attractive  for  studies  of  basic  circular  LGP  jet 
combustion.  This  is  because  the  pilot  shaft  protects  die  liquid  propellant  in  the  injector  from 
being  ignited  by  the  hot  ambient  gas  prior  to  injection.  Combustion  tests  with  LGP1845 
circular  jets  below  2  mm  in  diameter  have  been  conducted  recently^. 

4.  Thick  circular  jets  and  particularly  annular  jets  will  be  studied  primarily  in  non -combusting 
environment  using  high  speed  cinematography  and  flash  X-ray  radiography. 
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NOMENCLATURE 


-  Augment  action  area 

aelp  -  Exponential  arrbient  pressure  rise 

-  Hold  action  area 
A,  -  Inside  piston  area 

^  -  Effective  plugs  (part  28)  cross  section  area 

A^  -  Outside  piston  area 

^PS  -  Cross  section  area  of  pedestal  (part  25) 

Aj  -  Effective  test  chamber  gas  supply  inlet  area 

A^  -  Augment  volume  inlet  area 

Aj  -  Injection  port  area 

Ajj,  -  Hold  volume  discharge  area 

Cp  -  Discharge  coefficient 

M  -  Mass 

-  Piston  mass 

M  -  Molecular  weight 

P  -  Pressure 

R  -  Universal  gas  constant 

t  -  Time 

T  -  Tenperature 

Uj  -  Injection  velocity 

-  Piston  velocity 

V  -  Volume 

VL  -  Liquid  volume  injected 
y  -  Specific  heat  ratio 
PL  -  Liquid  density 

SUBSCRIPTS 

A  -  Atmospheric 

ALJ  -  Augment 

ADO  -  Augment  at  time  zero 

C  -  Test  chamber  (arrbient  conditions) 

H  -  Hold 

HO  -  Hold  at  time  zero 
HOP  -  Hold  opening  pressure 
L  -  Liquid 

S  -  Supply  reservoir 
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Intentionally  left  blank. 
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APPENDIX  A: 

INTERACTIVE  COMPUTER  CODE  INJFLOW 
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Intentionally  left  blank. 
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0  0,0 


APPENDIX  A:  INTERACTIVE  COMPUTER  CODE  INJFLOW 

program  INJFLOW 

integer  atc,athc,athx,atx,cdx,ii)inj,kJ±,ll,n,nmo,nouLrunge,wp 
parameter  (n=5,noutf6) 

logical  aa,ffp,ffpp,init,paps,plug,ppa,ppcm,ppns,ptrue,upno, 

1  uptrue.uptruel 

character  name*  12.part2*3,part4a*6.part3*  l,partla*  l,partlb*2, 

1  part5*  1  ,pait4b*  3,part2a*  1  ,pait3a*4 

double  precision  aamp,aaug,aexp,ah,ahmp,ai,aih,aiot/imstuK>,aoio, 

1  ar.as.asa.auath.athO.atth.atO.cd.co  1  ,co2,c  1  ,f(n),gam  .gamah, 

2  gamaug,gamp,gamt,gam  1  >gmo^moo,hJihIhhv^OJiTtmptmw,pa,pac, 

3  pacl,pag,paug,pc,ph,phi(n),phinn, phoop.pnop.pl, plo,psjho4ol, 

4  rollfsavey(n),tij,ts,uj,vaoo,vatth,vaug,vcno,vco,vn,xIy(n) 

variables  common  to  main  program  and  fen  subroutine 
common  /alpha/aamp,aexp,ah,aKmp,aoio,co  1  ,co2,c  1  ,gamah,gamaug, 

1  gamp.gamt.gam  1  ,gmo,gmoo,pa,pag,paps,paug,ppcm,ps)rol,upno, 

2  vaoo.vcho 


numeric  variables  common  to  main  program 
common  /beta/atth,atx,init,inj,pac,vatth,wp 


and  out  subroutine 


variables  common  to  main  program,  fen  and  out  subroutines 
common  /delta/aa,aiot,roll,uptruel 


aaug 

aexp 

A/UJ 

A 

"exp 

dvl 

f  (i) 

cfy(i) 

-~5T~ 

PS 

rho 

ah 

ah 

gam 

=  y 

ts 

ai 

= 

AI 

h 

=  At 

uj 

ams 

= 

ams 

np 

=  Mp 

vc 

ao 

= 

Ao 

mw 

=  M 

vh 

aps 

ss 

APS 

pa 

=  PA 

X 

as 

As 

paug 

-  p 

yd) 

asa 

= 

atau 

pc 

y(2) 

at 

s 

AI 

ph 

=  3, 

y(3) 

ath 

= 

ath 

phocp 

=  P 

y(4) 

cd 

= 

c 

'■'D 

Pi 

ii 

rW 

( 

y(5) 

•J 

vc 

:? 

=i 


Prompt  Area 
Liquid  Choice 

print*,’Please  choose  the  liquid  to  be  injected.’ 
print  101,  *  1)  LP  *  2)  Water  *,  ’  3)  Miscellaneous’ 
read*,kk 
if(kk.eq.l)  then 
rho  =  1.46 
partS  *  ’P’ 
elseif(kk.eq.2)  then 
rho  =  1. 
partS  =  ’W’ 
elseif(lck.eq.3)  then 
partS  =  M’ 

print30,  ’  Enter  the  density  of  the  liquid  (g/cc).’ 
read* /no 
else 
goto  4 
eixuf 


Chamber  Pressure  Mode 
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printVPlease  enter  chamber  pressure  mode.’ 
print  10,’  1)  Preset’,’  2)  Exponential  rise.’ 
read*Jck 

iffkk.eq.l)  ppcm  =  .false. 
if(kk.eq.2)  ppcm  =  .true. 
if(ppcm)  then 

print*, ’Please  enter  exponent’ 
read*,  aexp 
endif 

Injection  Mode 

print*, ’Please  enter  the  injection  mode.’ 

print  10,’  1)  Hold  discharge.’,’  2)  Constant  hold  pressure.’ 

read*Jdc 

iffkk.eq.2)  paps  =  .true, 
iffkk.eq.l)  paps  =  .false. 

Injection  Type 

print*  ,’Please  pick  the  number  corresponding  to  the  type  of  injec 
luon  to  be  modeled.’ 

print  101,*  1)  Circular  injection  -  Mode  K’,’  2)  Annular  injec 
lnon  -  Mode  L’,’  3)  Annular  injection  -  Mode  J’ 
farmat(lx,a/a) 
formatflx, a/a/a) 
read*4nj 

iffinj.eq.l)  ffp  =  .false. 
if(inj.eq.2j  fro  =  .tnie. 
if(inj.eq.3)  ffpp  =  .true. 


print  30,  ’Please  enter  estimated  value  for  CD.’ 

rcad*,cd 

cdx  =  cd*10. 

write(unit  =  part3,  fmt  =  '(il)’)  cdx 
if(cdx.eq,10)  part3  =  ’O’ 

Opening  Mode  For  Mode  K  Injection 

if(.not.(ffp.or.ffpp))  then 

print*.  Please  enter  the  mode  of  opening.’ 
print  101,’  1)  Pilot  shaft’,’  2)  No  pilot  shaft’,’  3)  Arb 
litrary  (only^if  hold  discharge) 

if(kk.eq.l)  then 
ppns  =  .false. 
part2a  =  ’S’ 
else 

ppns  =  .true. 
part2a  =  ’N’ 


if(kk.eq.3)  then 
if(paps)goto3 

pnnt*.  Please  enter  opening  hold  pressure  (psig).’ 
read*,  phoop 
ppa  =  .true. 
part2a  =  ’A’ 
endif 
endif 

Injection  Orifice  Diameter  For  Mode  K  Injection 
if(.noL(ffp.or.ffpp))  then 

print*,  Please  choose  the  size  of  the  injection  orifice  (inche 

Is).’ 

print  14,’  1)  .020’,’  2)  .040’,’  3)  .080’,’  4)  .125’ 

read*4tc 

if(atc.eq.l)  then 
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at  =  .02 

part4a  =  pait5//’.020’ 
elseif(atc.eq.2)  then 
at  =  .04 

part4a  =  part5//’.040’ 
elseif(atc.eq.3)  then 
at  =  .08 

part4a  =  part5//\080’ 
elseif(atc.e^.4)  then 

part4a  =  part5//\125’ 
else 
goto  7 
enau 
endif 

format  (lx,a/a/a/a) 
format  (lx,a/a/a/a/a) 

Opening  Mode  For  Mode  L  Injection 


if(ffp.or.ffpp)  then 
if(ffp)  then 

nrifirt  ’DIaapa  n< r 


print*, ’Please  pick  the  number  corresponding  to  the  mode  of 
1  injection  to  be  modeled.' 

print  101,’  1)  With  plugs.’,’  2)  Without  plugs.’,’  3) 
lArbitrary.’ 

read*,wp 
if(wp.ea.3)  then 

print*, Tlease  enter  opening  hold  pressure  (psig).’ 
read*,  phoop 
ppa  =  .true. 


Injection  Annular  Orifice  Diameter  For  Injection  Modes  J  and  L 

19  if(ffp)  print  30,  ’Please  alter  inner  annular  diameter  (between 
11.00  and  1.31  inches}.’ 

if(ffpp)  print  30,  ’Please  enter  outer  annular  diameter  (between 
1  1.00  and  1.308  inches).’ 
read  26,at 

if(aLlt.l.00.or.at.gL 1.308)  goto  19 
atx  =  at*  1000 
atx  =  atx- 1000 

write(part4b,  fmt  =  ’(i3.3)’)  atx 
iffwp.eq.l)  plug  =  .true. 
if(wp.eq.2)  plug  =  .false. 
part3a  =  W//part3  Mpart5//’ .  ’ 
iffplug)  part3a  =  ’P7/part3//part5//’.’ 
iftppa)  part3a  =  ’A7/part3//part5//’.’ 

atO  =  at 

Reservoir  Supply  Temperature 

print  30,  ’Please  enter  reservoir  supply  temperature  in  Celsius.’ 
read*,  ts 
ts  =  ts+273. 
pac  =  pa/14.7 

Constant  Hold 

if(pa^s)  then 

part2  =  ’HHH’ 
goto  11 
entfif 

Hold  Discharge  Orifice  Size 

9  _  jirint  30,  ’Please  choose  the  size  of  the  discharge  orifice  (inche 
print  8,  ’  1)  .020’,’  2)  1/32’,’  3)  1/16’,’  4)  wide  open’. 
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1  ’  5)  miscellaneous' 
read*,athc 
if(athc.eq.l)  then 
ath  =  .02 
part2  =  ’020' 
elseif(athc.eq.2)  then 
ath  =  .03125 
part2  =  ’031’ 
elseif(athc.eq.3)  then 
ath  =  .0655 
part2  =  '062' 
elseif(athc.eq.4)  then 
ath  =  .200 
part2  =  '200' 
elseif(athc.eq.5)  then 

print  30,  Please  input  discharge  orifice  size  (inches), 
read  26,ath 
athx  =  ath*  1000 
write (part2,  fint  =  ’(i3.3)’)  athx 
else 
goto  9 

25  enmf 
athO  =  ath 

11  at  =  atO 

ath  =  athO 
init  =  .true. 

c  Maximum  Test  Chamber  Pressure 

5  print  30,  ’Please  enter  chamber  pressure  in  psi.' 

read*,  ps 

c  Hold  Pressure 

print  30,  'Please  enter  hold  pressure  in  psi.’ 
read*,  ph 

c  Augment  Pressure 

if  (ppcm)  paug  =  0. 
if(.noLppcm)  then 

6  print  30,  'Please  enter  augment  pressure  in  psi.’ 
read*,  paug 

if(paug.ge.ps)  then 
print* 

print*. 'Caution,  augment  pressure  too  high.’ 
goto  6 
endif 
endif 

23  format(lx,a3) 

26  format(£5.3) 

30  format(lx,a) 

31  if(jiot.(ffp.ori^>p))  then 

lf^ppcm)  part  la  =  ’KE’ 
name  =  part  1  a//part2//part2a//part3//part4a 
else 

partlb  =  'LP' 
if(ppcm)  partlb  =  'LE' 
iftfnpp)  partlb  =  'JP' 
if(ffpp and. ppcm)  partlb  =  'JE' 
name  =  partlb//pah2//part3a//part4b 
endif 

c  Hold  Volume 


33  print*,’Please  enter  the  hold  volume  (cc).’ 

•  . .  -9.125',’ 


lii 


print  14,’  1)  231.25  (max)’,’  2)  129. 


3)  95.08’,’ 


*Jkk 
e 
vh 


if(kk.eq.l)then 
'  =  231." 
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4)  Arb 


30 


elseif(kk.eq.2)then 
vh  =  129.125 
elseif(kk.eq.3)then 
vh  =  95.08 
elseif(kk.eq.4)then 

print*,  please  enter  volume’ 
read*,vh 
else 

Joto  33 
f 

conversion  area 

data  pa,vaug,aaug,asa  .ah  ,gam,mw,h0,hl/1013000.0, 1.37, 16.67 ,0.0075 
139.1,L4,28.0,l.d-6,5.d-4/ 
pacl  =  14.7*pac 
ps  =  (ps+14.7)*pac 
ph  =  (ph+14.7)*pac 
paug  =(paug+14.7)*pac 
phop  =  (phoop+l4.7)*pac 

at  =  114159*(at/2.0)**2 
ath  =  ath*2.54 
ath  =  3.14159*(ath/2.0)**2 
if(.not.(ffp.or.ffpp))  then 
ao  =  40.53 
ai  =  17.65 
vco  =  1047.6 
mp  =  2750.0 
else 

iffffp)  at  =  8.68-at 
iflffpp)  at  *  at-5.067 
iftffo)  ao  =  45.60  -at 
iftffpp)  ao  =  40.53  -  at 
iftffp)  ai  =  22.72  -  at 
if(ffpp)  ai  =  17.65  -  at 
vco  =  829.72 
mp  =  3950.0 


viiuu 

pc  =  ps 
pi  =  pc 

more  variables  defined 

garni  =  (gam+1.0)/2.0 
gmo  =  gam/(gam-1.0) 
gmoo  *=  1.0/gmo 

cl  =  gam/(pn**(gmoo/2.0)*mw**0.5/156100.0)*ath 

col  =  sqrt(ps**gmoo*m  w/ts/83 140000.0) 

co2=  sqrt(pa«g**gmoo*mw)/156100. 

pag  =  pa**(jcim  1/gam) 

gamp  =  gam**0.5/gaml**(gaml/(gam-1.0)) 

gaml  s  1.0/gaml**gmo 

gmo  =  (2.0*gmo)**0.5 

aih  =  (ai-ao)/ah 

aah  =  -aaug/ah 

col  =  gam*  gmo/co  1  *as/aih 

co2  =  gam*gmo/co2*asa/aah 

gamt  =  (3 .0*gam- 1 ,0)/2.0/gam 

aiot  =  ai/at 

if(ffp)  aiot  =  (ai+at)/at 
atth  =  aiot*at/ah 
vatth  =  atth*vh 

rol  =  ai*rho/2.0/cd**2/mp*aiot**2 

roll  =  rol*mp/ai 

aoio  =  (ao-aO/mp 

ahmp  =  ah/mp 

aamp  =  aaug/mp 

vcho  =  (vco-aih*vh)/aih 


07,0.32, 
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c 


vaoo  =  (vaug-aah*  vh)/aah 
gamah  =  gam*ah 
gamaug  =  gam*aaug 

iffffp)  amsn  =  (ams*(ao+at)-(ai+at-aps)*(ai+at-ams))/ams/ah 
initial  values 


nmo  = 
uptrue : 
uptruel 
upno  = 


4 

=  .true. 

=  .false, 
.true. 


phin  =  0.0 
psin  =  0.0 
upin  =  0.0 
pninn  =  1 
u  =  0 
kk  =  0 
x  =  0.0 


call  out/x.yJcjiame) 
init  =  .false. 

40  if(aa)  goto  60 

k  =  rungefnmo,y,f,x,h) 
if(k.ne.l)  goto  50 

c  Calculations 


50 


call  fcn(x,y,f) 
goto  40 

ujupno)  f(5)  =  (y(5)-hh)/h 

if(uptruel.and.abs(f(l))/phinn.gt.2)  h  =  h/2. 
if(.not.paps)phinn  =  abs(ffl)) 
if(hhv-y(3).gt.y(3))  h  =  h/lO. 
hhv  =  y(3) 

if((.not.ffp.and.ppns).or.(ffp.and..noLplug))  then 
plo  =  y(l)+y(2)*aih+pacl*aah 
if(ppa)  plo  =  y(l)-phop 
zzv  =  {y(l)-pro)/pac-14.7 
iffppa)  zzv  =  phop/pac-14.7 
iffpaps)  zzzv  =  -plo/aah/pac-14.7 

if(.not.ffp.and..noLpnns)  then 


plo  =  yi 
zzv  =  y(2), 


.  ms 

25/pac-14.7 


zzv  =  yf  . 

iffpaps)  zzzv  =  -plo/aah/pac-14.7 
endif 
iffffp.and.plug)  then 
plo  =  y(l)-y(2)*amsh 
zzv  =  y(2)*amsh/pac-14.7 

iffpaps  .and..noLppcm)  zzzv  =  -(y(l)-y(2)*amsh)/aah/pac-14.7 
endif 
11  =  U+l 
ii  =  ii+1 
kk  =  kk+1 
if(y(l).ltpa)  then 
h=l.d-4 

iff  uptrue)  uptruel  =  plo.lt  le-20 
iffppcm)  goto  52 
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52 


34 


if((uptrae.and.uptrael).and.U.eq.l)  then 

print* ,’Not  enough  hold  pressure,  should  be  greater  than’,  zzv 
grint*,’  Please  reenter.’ 

Joto  11 
f 

if(paps.and..noLppa)  then 
plo  =  plo+paug*aah 
^lf^uptrue)  uptruel  =  plo.lt.le-20 

if((paps .and. uptrue)  .and.(.not. uptrue  1  .and. 11  .ge.2))  then 

print*, ’Not  enough  augment  pressure,  snoula  be  greater  than’, 
1  zzzv 

print*,’  Please  reenter.’ 
goto  11 
encuf 

if(uptrue.and.uptruel)  then 
h  =  hO 

print  34, ’INJECTOR  OPENS ’.’INJECTOR  OPENS’ 
uptrue  =  .not.  uptruel 

entT* 

format(lx/17x,al4,18x,al4/) 
if(abs(f(5)).gt.l.e6.and.uptruel)  then 
upno  =  .false, 
nmo  =  n 

j)rint*^NOW  INTEGRATING  5  EQUATIONS’ 


else 

upno  =  .true. 

nmo  =  n-1 
endif 

if(y(5).gt50.)  hi  =  l.d-4 

if(y(5).gt.500.)  hi  =  l.d-5 

if(ii.lLiO)  goto  40 

ptrue  =  abs(f(l)).gt.phin/l.l 

iffppcm)  ptrue  =  abs(f(2)).gLpsin/l.l 

imiptruel)  ptrue  =  abs(f(5)).ge.upin/l .  1  ,and.abs(f(5)).lt.  1  .e6 

iffptrue)  h  =  h*2. 

if(n.gt.hlj  h  =  hi 

pnin  =  abs(ffl)) 

upin  =  abs(f{5>) 

psin  =  abs(f(2)) 

u  =  0 


c  Output 

call  out(x,yJck,name) 
goto  40 
60  continue 

kk  =  kk+1 
call  out(x,yJck,name) 
print*  .’Injection  Time  =  ’,  x-tij 
print*,’Total  number  of  steps  =  ’,  kk 
end 


c  Runge-Kutta  Method  For  Solving  Differential  Equations 

integer  function  rungeln.y.fpcji) 
double  precision  y,£x,n  j)hi,savey 
integer  j,m,n 

dimension  phi(5),savey(5),y(5),f(5) 
data  m/0/ 
m  =  m+1 
goto  (l,2,3,4,5),m 

1  range  =  1 
return 

2  do  22  j  =  l,n 

{BS*! 

22  yO)  =  savey(j)+0.5*h*f(j) 
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x  =  x+0.5*h 
runge  =  1 
return 

do  33  j  =  l,n 
phi(j)  =  phi(jj+2.*f(j) 
y(j)  =  savey(j)+0.5*n*f(j) 
runge  =  1 
return 

do  44  j  =  1  ji 
pm  =  PhifefCj) 
yO)  =  saveyO)+h*fO) 

runge  -  1 
return 

do  55  j  =  l,n 

y(j)  =  saveyO'>Kphi(j>+f(j))*h/6. 

m  =  0 

runge  =  0 

return 

end 

Differential  Equations 

subroutine  fcn(x,y,f) 
integer  n 
parameter  (n=5) 

double  precision  x/(n),y(n),pth,plo/f ,vcc,gam  1  ,gmo,gmoo,c  1  ,pag, 

1  gamp,gamt,aiot,rol  jroll  ,aoio,ahmp,pa,ps  ,ah,co  l.vcho, 

2  gamah,cll,vaoo,co2,paug,gainaug,aamp,cccfaexp 
logical  aa,uptruel,upno,paps,ppcm 

common  /arpha/aamp,aexp,an,ahmp,aoio,col  ,co2,c  1  ,gamah,gamaug, 

1  gamp, gam  t,gam  1  ,gmo,gmoo,pa,pag,paps,paug,ppcm  ,ps  jol.upno, 

2  vaoo.vcho 

common  /delta/aa^iopoll.uptruel 

iffupno)  y(5)  =  sqrt((y0*aoio-y(l)*ahmp+y(4)*aamp)/rol) 


aa  =.true. 

y(5)  =  o.o 

goto  80 
entfif 

vcc  =  vcho  +y(3) 
iffpaps)  then 
f(l)  =  0.0 
else 

pth  =  y(l)*gaml 


pth  =  y(I)*gaml 

ff  =  gamah/y(3)*y(l)*y(5) 

ell  =  cl/y(3) 

•  _'V  A _ 


if(pth.gLpa)  then 

f(l)  =  -gamp*cll*y(l)**gamt+ff 

else 

if(y(l).lLpa)  y(l)  =  pa  _ 


UfyUl-it-paj  yfij  =  pa 
ccc  =  (y(l)/pa)**gmoo-1.0 
f(l)  =  -cl  l*gmo*y(l)**gmoo*pag*sqrt(ccc)+ff 
endif 
endif 

if(uptruel)  then 

u(.notppcm)  f(2)  =  col*y(2)**gamt/vcc’"((ps/y(2))**gmoo-l)** 

1  0.5+y(2)*gamah*y(5)/vcc 

ifrnpcm)  f(2)  =  aexp*ps*exp(-aexp*x) 

f(3)  =  -ah*y(5) 

vcc  =  vcc-vcho+vaoo 

f(4)  =  co2*y(4)**gamt/vcc*((paug/y(4))**gmoo-l)**0.5+y(4)*gama 
lug*y(5)/vcc 

f(5)  =  y(2)*aoio-rol*y(5)**2.0-y(l)*ahmp+y(4)*aamp 
else 

f<2)  =  0.0 

iffppcm)  f(2)  =  aexp*ps*exp(-aexp*x) 
f(3)  =  0.0 
f(4)  =  0.0 
f(5)  =  0.0 
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endif 
80  return 
end 


subroutine  out(x,yJck,name) 

integer  nouMiJck,114nj,wp,atx 

Seter  (nout=6,n=5) 

:  precision  y(n),x,pac,pl,uj,roll  ,aiot,atth,vatth,dvI 
logical  aa,init,uptruel 
character*  12  name 

common  /beta/atth,atx,init,inj,pac,vatth,wp 
common  /delta/aa,aiotjoll,uptniel 


if(aa)  goto  S00 

pi  =  y?2)+roll*y(5)**2 

uj  s=  y(5)*ak>t/100. 

II  =  Qdc+lVlO 


dvl  =  vattn-atth*y(3) 
if(init)  then 

wnte(nout,fmt=70) 


W1  IK^IIUUItlUU-  /  V/ 

1  ’LIO’  ’PIST’ 

*  pint  lM.’PRESS’.’PRESS’.'VOL’.’INJ’.’VEL’, ’PRESS’, ’PRESS’/V 


print*.’ - =====  .  . =  - 

1  s-  •  ■••  =  1  *  it  IBS' .  — 

print  100,ll^,y(4)^jac,y(l)/pac,y(3),dvl,uj,y(2)/pac,pl/p 
lac,V\5)/100» 

open  (2,file=name,recl=  10,form= ’  formatted  ) 
write(2, 105)ll,x,y(4)*  1  ,d-7,y(  1  )*  l.d-7,y(3),dvl,uj,y(2)*  1  .d-7, 
lpl*l.d-7,y(5)/l0O. 
endif 

if^n^lt^!ue^,y(4)/pac,y(l)/pac,y(3))dvl,uj,y(2)^ac,pl/pac,y 

1  write(2,105)Il,x,y(4)*l.d-7,y(l)*l.d-7,y(3),dvl,u|,y(2)*l.d-7, 

lpl*  1  .d-7,y(5)/100. 
endif 

100  format(lx,i4,lx47.5,lx,f7.2,lx48.3,3(lx,f7.3),2(lx,f8.2),lx,f6. 

13) 

105  format(i4,lx47.5,lx^(lxfn.3),3(lx,f7.3).2(lx.f8.3),lx/6.3) 

1 10  format(3x,al,4x,a4,5x,a4,3x,a4,4x,a4,6x,2(a3,4x),a5,6x,a3,4x,a4) 

120  format(16x,a5,2x,a5,5xIa3,6x12(a3,4x),a5,3x.a5,4x,a3) 
return 

500  close(2) 

501  return 
end 
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APPENDIX  B: 

FLOWCHART  OF  INERACTIVE  COMPUTER  CODE  INJFLOW 
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APPENDIX  B:  FLOWCHART  OF  INTERACTIVE  COMPUTER  CODE  INJFLOW 
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1  Commander 

US  Army  Materiel  Command 
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5001  Eisenhower  Avenue 
Alexandria,  VA  22333-0001 

1  Commander 

US  Army  Laboratory  Command 
ATTN:  AMSLC-DL 
Adelphi,  MD  20783-1145 

2  Commander 

US  Army,  ARDEC 

ATTN:  SMCAR-IMI-I 

Picatinny  Arsenal,  NJ  07806-5000 

2  Commander 

US  Army,  ARDEC 

ATTN:  SMCAR-TDC 

Picatinny  Arsenal,  NJ  07806-5000 

1  Director 
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and  Chemical  Command 
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1  Director 

US  Army  Aviation  Research 
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1  Commander  (Code  3247) 

Naval  Weapons  Center 
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US  Army  Materials  Technology 

Laboratory 
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1  Sandia  National  Laboratories 

ATTN:  R.  Rychnovsky, 

Division  8152 
P.O.  Box  969 

Livermore,  CA  94551-0969 

1  Sandia  National  Laboratories 

ATTN:  R.  Carling,  Division  8357 
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1  Sandia  National  Laboratories 

ATTN:  S.  Griffiths,  Division  8244 
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«  Physics  Laboratory 
fms  Hopkins  University 
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1  Jet  Propulsion  Laboratow 

ATTN:  D.  P.  Maynard,  Technical 
Group  Supervisor 
Chemical  Processes  Group 
Propulsion  Systems  Section 
4800  Oak  Grove  Drive 
Pasadena,  CA  91109-8099 

2  Princeton  Combustion  Research 

Laboratory 

ATTN:  N.  A.  Messina 
M.  Summerfield 
4275  U.S.  Highway  One  North 
Monmouth  Junction,  NJ  08852 

1  Calspan  Corporation 

ATTN:  Technical  Library 
P.O.  Box  400 
Buffalo,  NY  14225-0400 

7  General  Electric  Ordnance  System 

Division 

ATTN:  J.  Mandzy,  OP43-2202 
R.  E.  Mayer 
W.  Pasko 
R.  Pate 
I.  Magoon 
D.  Bair 
L.  Walter 
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General  Electric  Company 
Armament  Sytems  Department 
ATTN:  D.  Maher 
Burlington,  VT  05401 

Veritay  Technology,  Inc. 

ATTN:  E.  B.  Fisher 
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Pennsylvania  State  University 
Department  of  Mechanical 
Engineering 
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University  Park,  PA  16802 
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USER  EVALUATION  SHEET/CHANGE  OF  ADDRESS 


This  Laboratory  undertakes  a  continuing  effort  to  improve  the  quality  of  the  reports  it  publishes. 
Your  comments/answers  to  the  items/questions  below  will  aid  us  in  our  efforts. 

1.  BRL  Report  Number _ BRL-TR-3209 _ Date  of  Report  February  1991 

2.  Date  Report  Received _ _ _ _ 

3.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related  project,  or  other  area  of  interest 

for  which  the  report  will  be  used.) _ _ _ 


4.  Specifically,  how  is  the  report  being  used?  (Information  source,  design  data,  procedure,  source 
of  ideas,  etc.) _ _ _ 


S.  Has  the  information  in  this  report  led  to  any  quantitative  savings  as  far  as  man-hours  or  dollars 
saved,  operating  costs  avoided,  or  efficiencies  achieved,  etc?  If  so,  please  elaborate. _ 


6.  General  Comments.  What  do  you  think  should  be  changed  to  improve  future  reports?  (Indicate 
changes  tb  organization,  technical  content,  format,  etc.)  _ _ 
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City,  State,  Zip  Code 
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Address  in  Block  6  above  and  the  Old  or  Incorrect  address  below. 
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City,  State,  Zip  Code 
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